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Abstract 

In this paper we complete the results of Sullivant and Sturm- 
fels [SS] proving that many of the algebraic group-based models for 
Markov processes on trees are pseudo-toric. We also show in which 
cases these varieties are normal. This is done by the generalization 
of the discrete Fourier transform approach. In the next step, follow- 
ing Sullivant and Sturmfels, we describe a fast algorithm finding a 
polytope associated to these algebraic models. However in our case 
we apply the notions of sockets and networks extending the work of 
Buczyhska and Wisniewski [BW] who introduced it for the binary case 
of the group Z2. 

Introduction 

In the recent years phytogenetic trees have been paid much attention to 
because of the increased interest in evolutionary processes. They describe 
the DNA sequence changes and therefore they are a powerful tool that helps 
biologists to explain the evolution. They allow also for a better understanding 
of the relations between the species based on their DNA structure. The study 
of molecular phylogenetics is quite new but already advanced. 

Phylogenetic trees turned out to be very interesting also for mathemati- 
cians who recognized their relationship with algebraic varieties. Cavender 
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and Felsenstein |CFj as well as Lake [L] pioneered the work in algebraic phy- 
logenetics in the late 80s by introducing the invariants, that are polynomials 
describing a variety. 

In 1993 Evans and Speed |ES] observed that there is a natural group 
action of Z2 x Z2 on the nucleobases {A, C, G, T}, that defines models well- 
known from the biology. The method was generalized in |SSEj by proving 
that the obtained variety is pseudo-toricQ. Their results were based on the 
discrete Fourier transform and strongly relied on the hypothesis that the 
group acting on the space of states is abelian. Surprisingly they are valid 
also for some sub-models of the primary models. For example, the resulting 
varieties were known to be pseudo-toric for 2-Kimura and Jukes-Cantor mod- 
els even before the publication [SSJ. Still, it has not been specified exactly for 
which group-based models the variety is pseudo-toric. Further information 
on algebraic methods in phylogenetics can be found for example in [SSj or 
[PS] and the references therein. 

Many results in phylogenetic algebraic geometry, like the description of 
the generators in [SSj . were based on the pseudo-toricness blanket assumption 
on varieties. The arguments which imply that algebraic varieties related to 
submodels of general group-based models are pseudo-toric were commonly 
believed to work for any group-based model, although the authors agreed [Pr] 
that the assumption underlying their work may not be fulfilled and should 
have been stated more precisely. For an example in which the arguments 
implying pseudo-toricness of varieties related to submodels of general group- 
based models do not apply see the appendix. This problem seems to be 
crucial for the group-based Markov processes on trees that do not necessarily 
arise from the biology. 

Although the inspirations for phylogenetic models come from biology, 
they can be defined in a purely mathematical setting. We believe that, al- 
though mathematicians are often motivated by the empirical results, the 
main object of their study are still general theories. The fact that they are 
applied very often in other fields or in the real life is rather the special feature 
of mathematics than the aim of mathematicians. In case of phylogenetics the 
first models were introduced by biologists. Then a mathematical setting was 
developed in which biological models were covered. However other models 

In most papers concerning phylogenetics authors refer to pseudo-toric varieties as toric 
ones. In this paper we are going to discuses normality, so we will distinguish between these 
two notions. 
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(like a binary model) turned out to be also very interesting. The study 
of these objects resulted in many beautiful theorems and unexpected con- 
nections - see |SX] . |BW] . |Man] and |DK] . We believe that the setting of 
G-models introduced in this paper is a good foundation for developing the 
theory of group-based models. They are sufficiently general to cover all the 
group-based models of interest and still, as it has been shown, have all the 
expected properties. 

The main aim of this paper is to prove that a large class (including well 
known 2-Kimura and Jukes Cantor models) of group-based algebraic models 
are pseudo-toric. We will do this by using not only abelian groups, but also 
arbitrary groups that have a normal, abelian subgroup. 

The main theorem of this paper 13.191 proves the pseudo-toric result in a 
general setting: 

Main Theorem. Let H be a normal, abelian subgroup of a group G C S n = 
Sym(A), that acts transitively and freely on a set A of n elements. Let W be 
the space of matrices invariant with respect to the action of G and let W be 
the vector space spanned freely by elements of A. Then the algebraic model 
for (T, W, W) is pseudo-toric for any tree T. 

Moreover in such a setting the main theorem 26 of |SS] is applicable. 

Furthermore we investigate which varieties are normal and thereby toric. 
For necessary information concerning toric varieties see |Ful] . |Oda] or |CLSj . 
It turns out that many examples arising from abelian groups (like 3-Kimura 
model) are toric, but for arbitrary groups (the case of 2-Kimura) this is not 
true, even in the case of the simplest 3- leaf tree. 

This paper does not demand any knowledge from nonmathematical sci- 
ences. Readers interested solely in algebraic geometry can skip remarks con- 
cerning biological or probabilistic setting, although the intuition form Markov 
processes can be helpful. 

The structure of the article looks as follows. In the first section we recall 
basic facts about constructing varieties associated to the phylogenetic trees. 
We will try not only to state precise definitions but also to give some intu- 
itive meaning of introduced objects. In the second section we present results 
concerning group-based models in the abelian case. We also restate and con- 
ceptualize notions of sockets and networks from |BWj . In section three we 
define G-models that are generalizations of the abelian case. In a slightly 
different setting the idea of G- models can be found in |BDWj . The family of 
(7-models contains such examples as 2-Kimura, Jukes-Cantor and 3-Kimura 
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model. Our results apply not only to biological phylogenetic models but also 
to many other Markov processes. In section four we present efficient algo- 
rithms used for constructing polytopes of the pseudo-toric models. In section 
five we prove facts concerning normality of a model and we show examples 
in which the obtained variety is not normal. To prove normality we reduce 
any tree to the simplest 3-leaf tree checking it by the use of the Macaulay 
computer program. In the last section we state some open questions. 
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1 Preliminaries 

In this section we will introduce a general setup for the rest of the paper. We 
remind basic facts about varieties associated to phylogenetic trees. From the 
point of view of Markov processes each point of our variety can be considered 
as a generalization of a probability distribution on the space of states of 
leaves of the tree. First we fix a rooted tree T. In many biological cases 
this tree encodes the history of the evolution of DNA. In our setting we 
assume that the root distribution is uniform. The root is thought of as a 
common ancestor, each edge corresponds to mutation and leaves correspond 
to species that mutated the last, and so are observable. If we speak about 
directed edges of a tree T we always assume that they are directed from the 
root. However it turns out that the variety associated to a given tree does not 
depend on the orientation of edges (so neither on the root) for any (7-model, 
see remark 12.41 

We fix a finite set A. In the biological setting the set A corresponds to 
the possible observed states on a given position in the DNA sequence. In 
biology the elements of A are denoted by A, C, G, T. We will be considering 
more general Markov processes, so A can be an arbitrary finite set. 
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1.1 Algebraic setting of phylogenetic trees 

In biology we assume that we do not know the state of inner vertices of 
trees. That is why we associate to each node an A-valued random variable. 
Of course such variables give as a distribution - a set of nonnegative real 
numbers that sum up to 1. Following |ERSSj we generalize this setting. 
We work over the field of complex numbers C and we omit the condition 
of summing to one. In this setting in each vertex v of the tree we obtain a 
complex linear space W v = W with the distinguished basis given by elements 
of A. Random variables described above correspond to the set of points with 
real, nonnegative entries that sum up to one - this is a simplex. 

Apart from the states, for a specific model we also define the possible 
transition mechanism along edges (like a Jukes-Cantor model or Kimura 
models). In biology this part of the data depends on the model of the muta- 
tion mechanism that we choose. In general, for a Markov process it is defined 
by an association of a bi-stochastic matrix to each edge. The entries of this 
matrix correspond to relative probabilities of mutations along edges. In an 
algebraic setting (see also |ERSSj ) the bi-stochastic condition is replaced by 
the condition that the sum of all coefficients in columns and rows are equal. 
For now, let us forget about this condition and we will see later that in our 
case it will be always satisfied. For the time being we associate to each edge 
from v\ to V2 a matrix, that is a linear map from W V1 to W V2 in the chosen 
basis. As in each vertex we have a distinguished basis, we can view the space 
of morphisms form W Vl to W V2 as W Vl g) W V2 , namely: 

End(W Vl ,W V2 ) 3 (f : W V1 W V2 ) -> ^ <7 2 *(/(<?i))<7i ® 92, 

(91,92)&G V1 xG V2 

where G v . is the distinguished basis of W Vi and g* is an element of the basis 
dual to In a given phylogenetic model we consider only matrices from 
some subspace W of W <g> W. Often we assume for example that the matrix 
is symmetric or as already mentioned bi-stochastic in an algebraic sense. 

The rest of this section fixes some notation and can be found in |BW] . 
NowJ^o each edge e of a given rooted tree T we associate a copy of W denoted 
by W e (in fact, we can associate different vector spaces to different edges and 
all the arguments will work). 

Definition 1.1. We define the space of all possible states of a tree: 

W V = ®W V , 

vev 
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where V is the set of vertices of a tree, the parameter space: 

W E = (g)W e , 

eG-E 

where E is the set of edges and the space of the states of leaves: 

W L = (£)W h 

leL 

where L is the set of all leaves of the model. 

Let (vi,v 2 ) be an edge that connects two adjacent vertices v± and v 2 
(directed from v\ to v 2 ). We consider a linear map: 

(W V1 <g> W V2 ) -»• W v . 

Intuitively this map associates to a particular choice of matrices over all 
edges the "probability distribution" on the set of all possible states of the 
tree. In order to define this map we have to define it on the basis. As A can 
be considered clS 8b basis of W the basis of 

(g) (W V1 ®W V2 ) 

(v 1 ,v 2 )€E 

can be chosen as follows: 

e=(vi,v 2 )£E 

where for each e adjacent to v the elements a v e are elements of basis of W v . 
Such an element corresponds to setting a matrix over an edge e = (vi,v 2 ) 
that gives all probabilities of mutation equal to apart from the probability 
of changing from a^ 1 to a^ 2 along edge e equal to 1. Now it is easy to set the 
corresponding probability on the set of states: 

?(®e=(« 1) « a )es(ar ® a e 2 )) = ®veva v , 

if for any v all a v e are equal for all edges adjacent to v (we define a v as the 
vector that they are all equal to) and otherwise. 
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Definition 1.2. We define the map: 

$: W E -> WV, 

as the restriction of if)' . This corresponds to the restriction of the possible 
space of matrices over any edge. 

The closure of the induced rational map 

iP:l[F(W e ) ~+F(W v ), 

is called the complete projective geometric model. 

Inspired by biological setting mathematicians consider not the complete 
projective geometric model, but a model that takes into account that we do 
not know the states of inner vertices of the tree T. We define a map: 

5 : W -> C, 

such that 5(a) = 1. In other words 5 = J2 a &A a * ■ Next we define a map that 
hides the inner vertices. 

Definition 1.3. Let 

vr L : W Y -»■ W L , 

be a map defined as 

where the map 5w v acts like 6 on the vector space associated to a node v. 

Remark 1.4. The idea of the map ttl is the following: Consider points of 
Wy as generalizations of probability distributions on the space of all possi- 
ble states of the tree T and points of Wl as generalizations of probability 
distributions on the space of possible states of the leaves of T. Let p e Wy. 
The image 7r(p) corresponds to the distribution that to a given state S of 
leaves associates the sum of all probabilities (induced by p) of states of the 
tree T that coincide with S on leaves. This remark can be read literally if 
we restrict the spaces Wy and Wl to simplices described at the beginning of 
the section, where the points have got a probabilistic meaning. 

If we compose the map ip with ttl we obtain a map from We to Wl- 
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Definition 1.5. The map described above induces a rational map: 



\F(W e ) -+F(W L ). 



The projective geometrical model, or just a model, denoted by X(T) is the 
closure of the image of this map. 

2 Group models - abelian case 

Most of this section is well known. However we repeat some of the argu- 
ments in order to fix the notation and introduce some ideas that will be 
needed in following sections. We also generalize the notions of "sockets" and 
"networks" introduced in |BWj . what enables us to generalize some of the 
results form Z2 to arbitrary abelian groups. 

Let T be a rooted, trivalent tree. Let A be a set of letters (possible states). 
In their famous paper [ESJ Evans and Speed recognized a natural action of 
an abelian group G on A in biological case, namely the group G = Z2 x Z2 
acts on {A, C, G, T} transitively and freely. Their method was described in 
a more general setting in |SSEj . From now on we assume that we have a 
transitive and free action of an abelian group G on A. Of course, if we define 
the vector space W with the basis A, then the action of G on A extends to 
the action of G on W. 

Definition 2.1. For g G G let A g be a matrix (equivalently a linear map) 
corresponding to the action of g on W. 

By choosing one element of the set A and associating it to a neutral ele- 
ment of G we may make an action preserving bijection between the elements 
of A and G. An element associated to a G A will be denoted by g a . This 
allows us to find another basis of W, indexed by characters of G. This is 
done by the discrete Fourier transform. 

Definition 2.2. Let x £ G* be any character of the group G. We define a 
vector w x EW by: 




ae A 
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One can prove that the elements w x form a basis of W. Let us notice 
that although the choice of the bijection between A and H is not cannonical, 
the one dimensional spaces spanned by w x are. Changing the bijection just 
multiplies each vector w x by xio) f° r some g e G. 

The group structure distinguishes also naturally a specific model, namely 
the vector space W. The space W C End W, but as we have a distinguished 
basis on W made of the elements of A, we have an isomorphism W* ~ W, so 
as already mentioned we may identify End W ~ W ® W. We have a natural 
action of G on W Cg> W - the action of g is just g ® 

Aai g) a 2 ) = Ay(ai) <g> #(a 2 ). 

Now, we may define W as the set of fixed points of this action - see also 
jBDW] . 

Remark 2.3. In other words we take only such matrices that satisfy the 
following condition for any g G G: 

If we permute the columns and rows of matrices with a permutation cor- 
responding to g then we obtain the same matrix. 

Remark 2.4. One can see that if A e W, than A T e W. This means that 
if we consider a tree T with two different orientations than the models are 
exactly the same. If a point is the image of some element of the parameter 
space with respect to a given orientation than it is also the image of an 
element of the parameter space with respect to the second orientation. We 
just have to transpose matrices that are associated to edges with different 
orientation. 

Example 2.5. Let us consider the group Z 2 = (0, 1). We obtain matrices of 
the form: 

a b 
b a 

This is a binary model well studied in IBWf . 

Example 2.6. In this example we consider the case of particular interest in 
biology. The natural action on A = {A, C, G, T} is the action of the group 
Z 2 x Z 2 . Let us consider the group G = Z 2 x Z 2 . We obtain a model: 



a 


b 


c 


d 


b 


a 


d 


c 


c 


d 


a 


b 


d 


c 


b 


a 
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that is a 3-Kimura model. 



One can also check that using the basis (w x ) e Q the elements of W are 

exactly the diagonal maps. This means that we can define elements of W 
also indexed by characters of the group G: 



Definition 2.7. 

l x (w x > 



w x X = x' 

x + x' 



It follows that these elements are a basis of W. One can also check that 
using the basis A we have: 

Proposition 2.8. 

We see that the vectors l x are independent from the choice of the bijection 
between A and G, as the element g^}g a is a unique element of G that sends 
a to a, hence does not depend on the bijection. The map l x is a projection 
onto the (canonical) one dimensional subspace spanned by w x . 

Using this basis we will see that the map defining the model is given by a 
subsystem of the Segre system. First we need some lemmas. Of course, the 
action of G on W extends to the action of G on Wy and Wl- We have: 



Lemma 2.9. 



di^Wy 4 = ^\G\\ y \ = \G\\ y \- 1 



\G\ 



&mW? = ^\G\\ L \ = \G\\ L \~ l 
Proof. Let us consider the following basis of Wy- 

The action of g in this basis is diagonal, so the space of invariant vectors is 
spanned by invariant elements of this basis. As g{w x ) = x{g~ l ) w x we obtain: 

g(®vevW Xv ) = ®ve_vXv{g~ l )w Xv = Yl Xvig' 1 ) ®vev w Xv , 
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so such an element is invariant iff for any g £ G we have Husw^O?) = 1- 
This is equivalent to the condition that ^2 veV Xv is equal to a trivial character 
(we use additive notation for the group G*). From this we see that the 
dimension dim Wy is equal to the number of sequences indexed by vertices 
of the tree of characters that sum up to a neutral character. This gives 
us sequences what proves the first equality, as for abelian groups 

|G| = \G\. The proof of the second equality is the same. □ 

Remark 2.10. The basis {®„ 6 yw x „} of Wy in not natural (that is depends 
on the choice of the bijection between the set A and G). However the basis 
{® V £V w Xv : J2v€V Xv = Xo} of Wy is natural. Changing the bijection mul- 
tiplies w x by x(q) f° r a fixed g £ G, so as J2vev Xv = Xcb then of course 
(E veV Xv){a) — 1 and the vectors remain unchanged. 

One can easily see that the image of We in Wy is invariant with respect 
to the action of G. 

Proposition 2.11. The following morphism: 

ij:HF(W e )^F(WC), 

eG-B 

is given by a full Segre system. 

Proof. We prove that there is an isomorphism of We and Wy. Thanks to 
lemma ESS the dimensions are the same, because the number of edges is equal 
to number of vertices minus one, so it is enough to prove the surjectivity of 
the dual morphism W v — > W E . In this proof we always consider the basis of 
W given by elements of A. The space W E is generated by tensor products 
of functions that for a given matrix in W return a given entry of the matrix: 

® eeE (a' e )* (8) «)* 

where a* is an element dual to a and if e = (vi,v 2 ), then a' e £ W V1 and 
a" £ W V2 are elements of the chosen basis. Let g e be a unique element of G 
that sends a' e to a". Next, we construct an element of Wy inductively starting 
from the root, where we take any element from A. Suppose that we have 
already defined an element a for a vertex v. We take an edge e = (v, w) and 
we define the state at w as g e a. One can easily check that the image of the 
dual of the tensor product of so defined elements gives us <X>eG-E( a e )* ® ( a e)*> 
what proves the theorem. □ 
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Moreover the isomorphism of We and Wy has got a very nice description 
in the new basis. 

Proposition 2.12. The isomorphism of vector spaces We and Wy takes 
the base {® eeE \G\l Xe } bijectively onto the base : ^2 veV Xv = Xo}, 

where Xo is a trivial character. This means that in these basis the Segre 
embedding is given by monomials. 

Proof. We can see that: 

1 1 . 

' ' ' ' aGA 

This shows that the image of ® eeE \G\l Xe has got c = Ue^^&E Xe(9a^9a V2 ) 
as a coefficient of ® v& ya v . The coefficient c is equal to: 

c= JJ (-Xe)(g avi )xe(g av2 )- 

For given characters Xe let us define characters Xv f° r & U v vertices of the 
tree as: 

Xv = 

(v,w)£E (w,v)eE 

This corresponds to summing all characters on edges adjacent to v with 
appropriate signs, depending on the orientation of the edge. Using this no- 
tation: 

vev 

We define an element ® v evw Xv that is clearly in the chosen basis of Wy 
as each character Xe is taken twice with different signs, so the sum of all Xv 
is a trivial character. Moreover 



aGA 

so the coefficient of (^v^vo-v is 

equal to Yl veV Xv(9a v ), what proves the theo- 
rem. □ 

Let us notice that apart from the action of G on W ® W given by g ® g 
that allowed us to define W, we have got another action of G on W ® W 
given by g <S> id, where id is the identity map. 
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Lemma 2.13. The action g ® id restricts to W. 

Proof. It is enough to prove that the image of the action of g ® id on any 
element that is invariant with respect to the action g' <g> g' is also invariant. 
Let G be any element of W. 

(</ ® </)((<? ® id)c = (</<? ® </)(<?) = (<?</ ® = (g ® «W ® </)(<?) = 

(<7<gud)(G). 

Here we used the fact that G is abelian. □ 

Definition 2.14. Let N be the set of nodes of the tree. We define fP v for v 
an inner vertex and g G G as an isomorphism of the space We that acts on 
edges (v, w) as g <8> id, on edges (w, v) as g~ x ® id and on edges not adjacent 
to v as identity. We also define a group Gn as a group generated by ffl v . 

Remark 2.15. It is crucial to realize how g <g> id acts on elements of W 
considered as morphisms. One can check that g <g> id(A g /) = A g > o A g -i, so 
the action of g £g> id composes the given morphism with A g -i. 

Remark 2.16. The definition of p% corresponds to changing the matrices 
above edges adjacent to v in such a way that they give the "probability" as 
if the state of v was different. We will see in the following proposition that 
these definitions are useful, because if we do not know the states of inner 
vertices, then the orbits of the action of Gat give us the same "probability 
distribution" on the space of states of leaves. 

First we need a technical lemma. 

Lemma 2.17. The group G N ^ GM. G n acts on W E . There is a base in 
which Gn acts diagonally. 

Proof. We will prove that the good base is the tensor product of endomor- 
phisms l x . Using we obtain: 

(g ® id{l x )){w x >) = l x A g -i(w xl ) = 

= ^(E^W ) = l xC^Z^(ga)g~ l o) = 
= l xC£2x'(g a g)a) = x'(g)l x ( w x') = x(g)i x ( w x')i 

where the last equality follows from the fact that l x {w x >) is non zero only if 
X = x' ■ This proves that g ® id(l x ) = x(g)Z x , what proves the theorem. □ 
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Remark 2.18. As described above the elements of the base of We are bi- 
jective with the sequences of characters indexed by edges of a tree. In other 
words an element of a basis of We can be described as an association of a 
character of G to each edge of a tree. Moreover the elements of the basis 
of We that are invariant with respect to the action of Gn are exactly such 
associations that a (signed) product of characters around each inner vertex 
gives a trivial character. Precisely, for each vertex v and edges eo = (wo,v), 
ei = (v, Wi), . . . , &i = (v, Wi), . . . , e k = (v, w k ) we have: 

k 

~Xeo ^ ] Xe k — XOj 
i=l 

where the sum is the sum in G (remember about the additive notation in 
G), Xe is the character associated to an edge e and %o is a trivial character. 

Proposition 2.19. The map: 

W v 3 ® v€V 9v -> ®ieL9i e W L 
in the basies w x is given by: 

w v 3 ® vev w Xv i^r 1 ®ie£ w xi 

if all the characters for the inner vertices are trivial or zero otherwise. 
Proof. First let us look at ® v e_v w \v i n the old coordinates: 

where the sum S(o»)„ yeA v * s taken over all | V|-ples (indexed by vertices) of 
basis vectors. In other words this sum parameterizes the basis of Wy made 
of tensor products of base vectors corresponding to elements of G. This is of 
course equal to: 

Y Y II Xv (ga v ) JJ Xf (9a f ) ®v€N d v ® /eL df . 

(a u ) ueN £A N ( ai ) leL £A L v£N f<=L 

We see that the image in Wl is equal to: 

Y Y II Xv(g av ) JJ Xf(g af ) ®feL a f = 

(a u ) ueN £A N (at) leL eAL v€N feL 
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The product (Y[ u eN(^2geG Xu(g))) is equal to zero unless all characters Xu 
for u E N are trivial. In the latter case the product is equal to iGj'^L Of 
course 

(II Xf(gi))(®ieLgi) = ®leLW xv 
(gi)ieL£G" feL 

what proves the proposition. 

□ 

The following theorem is a direct generalization to arbitrary abelian 
groups of theorem 2.12 from |BWj . 

Theorem 2.20. The spaces (W G ) and (We) Gn are isomorphic. 

Proof. One can prove it using dimension argument, but it is better to look 
how the basis are transformed. The base of (We) Gn is given by <S) e eE\G\l Xe , 
where the signed sum of all characters at any vertex is trivial. This, thanks 
to I2.12[ by the morphism We Wv is transformed (bijectively) into an 
independent set ® v evW Xv , where characters for inner vertices are trivial and 
the sum of all characters is trivial. Using 12.191 the image of this set via 
Wy — > Wl gives us (bijectivly) \G\^ ®«eL w xv where the characters sum up 
to a trivial character. The last set is the basis of W/J. □ 

This motivates the following definitions of sockets and networks. 

Definition 2.21. A socket will be an association of characters from G to 
each leaf such that the sum of all these characters is a trivial character. 

A network will be an association of characters from G to each edge such 
that the (signed) sum of characters at each inner vertex gives the trivial 
character. 

Let us generalize the results on sockets and networks from |BWj . One 
can see that there is a natural bijection of networks and sockets: a network is 
projected to leaves and a socket is extended to a network using the summing 
condition around inner vertices. Moreover each network determines naturally 
an element of the basis of (We) Gn and each socket an element of the basis 
of WfJ. The isomorphism of 12.201 just uses this natural bijection. Of course 
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as G was abelian instead of associating characters we can associate group 
elements. 

In case of the group Z 2 ( |B W] ) sockets were even subsets of leaves. That 
was associating 1 to chosen leaves and to the other leaves. The condition 
that the subset has got even number of elements is just the condition that 
the elements from the group sum up to a neutral element. We see that this 
definition is compatible. Networks were subsets of edges such that there was 
an even number chosen around each inner vertex - this is also a condition of 
summing up to a neutral element around each inner vertex. 

Using the theorem 12.201 we know that the model X is the closure of the 
image of the rational map: 

£\G\\E\ _^ p(|L|-l)x|G| 

where the coordinates of the domain are indexed by pairs (e,x) for e an edge 
and x G G*. The coordinates of the codomain are indexed by sockets (or 
equivalently networks). Let M be the lattice with the basis given by pairs 
(e, x)e&E, x &G*- Let P be a polytope that is a subpolytope of a unit cube and 
whose points correspond to networks. It follows that the model is given by 
the polytope P. 

To each pair (e, %) where e is an edge and x is a character of G we can 
associate a one parameter subgroup that is given as a morphism from M to 
Z and is a dual vector to the vector of the base of M that is indexed by the 
pair (e, x)- 111 particular for each leaf I and character x G G* we obtain a one 
parameter subgroup Xf. For each t G C* the action of Xf(t) on the model 
extends to the action on pd L l -1 ) x l G l 3 X. The weight of this action on the 
coordinate indexed by a socket s is either or 1 depending on whether the 
socket s associates to the leaf I character x (i n this case 1) or not (in this 
case 0). 

Remark 2.22. In |BW] the authors considered only one one parameter sub- 
group for each leaf although their group had two elements. Notice however 
that in our notation for the group Z 2 the weights of the action of A° are 
completely determined by the weights of the action of A; 1 - one weights are 
negations of the others. In our notation the authors considered only A; 1 . 
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3 Nonabelian Case 



The setting of this section is sufficiently general to cover many Markov pro- 
cesses, in particular this will be a generalization of the results of the previous 
section. However the inspiration is a 2-Kimura model, that is a phylogenetic 
model in which the transition matrices are of the following type: 



In this case, as in the previous section, we also have an abelian group 
H = Z 2 x Z 2 that acts on the basis (A,C,G,T) of a four dimensional vector 
space W. As we have seen the fixed points of the action of H on W <g> W 
define a 3-Kimura model. We may however define a larger group G, namely 
.Dg that contains if as a normal subgroup and the action of G on W ® W 
defines the 2-Kimura model (compare with [BDWJ). That is why we consider 
the following setting. 

Let A be an ro-element set of letters. Let G be a subgroup of S n = Sym(A) 
(not necessarily abelian) acting on A. Suppose moreover that the group G 
contains a normal, abelian subgroup H and the action of H on A is transitive 
and free. Elements of A once again correspond to states of vertices of a 
phylogenetic tree T. We define If to be a complex vector space spanned 
freely by elements of A, namely W = (BaeA^a, where C a is a field of complex 
numbers corresponding to one dimensional vector space spanned by a G A. 

What changes is the definition of W. We define elements of W as matrices 
fixed not only by the action of H, but by the whole action of G. We use the 
notation assuming that End{W) C W <g> W. 

Definition 3.1. Let 



Remark 3.2. The characterization of W from 12.31 is still valid. 

Remark 3.3. The situation of the previous section corresponds to G — H. 



a b c b 
babe 
c b a b 
b c b a 
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Remark 3.4. As before by choosing some element of A we may make a 
bijection between A and H. However this time we have to be very careful. 
The action of G on A (as permutation) will not generally be the same as 
the action of G on H (as a group). We fix one such bijection. An element 
associated to a £ A will be denoted by h a £ H. 

We will often use the following easy observation: 

Lemma 3.5. Let h £ H be such an element that as a permutation sends a 
to b, where a,b £ A. Then h = hbh~ l . 

Proof. Of course both elements send a to b, so because H acts on A transi- 
tively and freely, they have to be equal. □ 

Definition 3.6. A G -model, will be a phylogenetic model associated to a 
phylogenetic tree (T, W, W) . 

Our aim is to prove that also in this generalized setting we will obtain 
pseudo-toric varieties. We will proceed in four steps. 

1. We introduce a general method for constructing endomorphisms of W 
from complex functions on H. We prove that under certain conditions 
(namely a function should be constant on orbits of the conjugation 
action of G on H) the obtained endomorphism is in W. 

2. We prove that some sums (over the orbits of the action of G on H) of 
characters of H are functions that can define elements of W. We also 
notice that we obtain a set of independent vectors of W. 

3. Using dimension arguments we prove that the set defined in step 2 is 
in fact a basis. 

4. Finally, using theorems from section 3, we prove, using the new coor- 
dinates that our variety is pseudo-toric. 

Remark 3.7. Of course we can define Wh as a vector space of matrices fixed 
by the action of H. From the previous section we know that the closure of 
the image of the map: 
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is a pseudo-toric variety. Moreover we also found the base in which the 
described morphism is given by monomials. Of course W C Wjj, so our aim 
is to prove that the restriction of the previous map is also given by monomials 
in certain base. We will use the base on Wh to define the base of W. 

3.1 Correspondence between functions on H and en- 
domorphisms of W - step 1 

We are going to define some endomorphisms of W. 
Definition 3.8. Let f : H — >■ C be any function. We define: 

h=J2 f(K l h b )a®b. 

a,beA 

Remark 3.9. Notice that because of I2.8l the previous definition is consistent 
with the definition of l x for x £ H. Moreover the vector If depends only on 
the function / and not the bijection between A and H, as h~ l h b is the only 
element from H the sends a to b. 

Proposition 3.10. Let us consider an action of G on H: 

(g,h) ->■ ghg~\ 
If f is constant on orbits of this action then If G W. 

Proof. Consider any element g G G. We focus on two entries of the matrix 
If, namely (ai,&i) and (a2,b 2 ), where 

g{ai) = a 2 and g{b x ) = b 2 . 

These entries are from the definition of If respectively fih^h^) and f(h~^hb 2 ) 
so because of remark EOl we want to prove that: 

f(h^h bl ) = f(h~ 2 1 h h2 ). 

Consider an element gh bl h~^g~ l . Clearly it is an element of H (because H 
was a normal subgroup of G) that sends a 2 to b 2 . From lemma EJ31 we obtain: 

Because / was constant on orbits of the action of G and H was abelian we 
get: 

f(K'h bl ) = f(h~ 2 1 h b2 ), 
what completes the proof. □ 
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3.2 Appropriate functions on H - step 2. 

In the abelian case we considered characters of H. As G was equal to H, 
these functions were of course constant on (one element) orbits of the action 
of G on H. In a general case it may happen that we do not have an equality 

xighg' 1 ) = x(h)- 

Of course this equality holds if a character of H extends to a character of 
G, but this is not always the case. If we define the vectors l x for \ £ H 
they may not be in W . To obtain the vectors in W we will sum up some 
characters to obtain functions that satisfy the condition of 13.101 Consider 
the action of G on i/: 

(g,x)(h) = xighg' 1 )- 

Let O be the set of orbits of this action. Of course elements of O give a 
partition of H. Let us define for each element o £ O a function f : H — > C. 

Definition 3.11. Let 

Here we are summing characters as complex valued functions, not as char- 
acters, so this is a normal sum, not a product. We obtain 

XGo 

Proposition 3.12. The function f satisfies the conditions of \3.1U\ (is con- 
stant on orbits of the action of G on H by conjugation). 

Proof. 

Ug'hg'- l ) = C£x{g'hg'~ 1 ) = 

X6o 

= J2(9',x)(h) = J2x(h) = fo(h), 

X£o xSo 

because the action of g' is a permutation of the orbit o. □ 

Corollary 3.13. The vectors lf a for o £ O are in W. Moreover, as l x 
formed a basis ofWn, and lf o are sums over a partition of this basis, they 
are independent. 
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Proposition 3.14. Any function f that is constant on orbits of O is a linear 
combination of functions f Q . 

Proof. Of course / can be uniquely decomposed into sum of characters of H: 

We have to prove that coefficients of x m the same orbit are the same. Let 
(flS Xi) = X2- We know that for any h € H we have 

Y^a^xih) = f(h) = fighg' 1 ) = 

Yl a xx(ghg~ l ) = ^ a x(9, x)(h). 

From the orthogonality of characters we see that a Xl = a X2 what completes 
the proof. □ 

Corollary 3.15. The number of orbits in O (and so the number of vectors 
lj o ) is equal to the number of orbits of the action of G on H: 

(g,h) ->■ ghg~ l . 

Proof. This follows from comparing dimensions of spaces of complex func- 
tions on H that are constant on orbits. □ 

3.3 Dimension of W - step 3 

We are going to prove that the dimension of W is equal to the number 
of orbits \0\. First let us note that all coefficients of any matrix in W 
(in the basis A) are determined by coefficients in the first row (this follows 
from the second section). We see that dimly is equal to the number of 
independent parameters in the first row. Let e be a fixed element of the set 
A (corresponding to the first row of the matrix). The action of G imposes 
some conditions, namely the coefficient in the e-th row and a-th colomn and 
the coefficient in the e-th row and 6-th colomn for a,b G A have to be equal 
iff there exists an element g e G such that: 

g(e) = e and g(a) = b. 

By h c we will denote the (unique) element of H that sends e to c. 
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Lemma 3.16. The following conditions are equivalent: 

1. 3 9e G that sends e to e and a to b, 

2. the elements h a and hi, are in the same orbit with respect to the action 
(g,h) = ghg-\ 

Proof. Of course h a and hf, are in the same orbit iff h~ l and h^ 1 are in the 
same orbit. For the proof we concentrate on the second variant. 

1 =>- 2: From 1X51 we know that gh~ 1 g~ 1 = h^ 1 , because both elements 
send b to e. 

1 <= 2: Suppose that gh~ l g~ l = . Let g' = gh g -iQ>). Of course g' 
sends e to e, but g' = gh^hg-ia,), hence it also sends a to b. □ 

Proposition 3.17. The dimension of W is equal to the number of orbits 
\0\. 

Proof. Equal parameters in the first row of matrices in W correspond bi- 
jectively to orbits of the action of G on H from 13.161 and remarks at the 
beginning of this subsection. This, along with I3.15[ finishes the proof. □ 

Corollary 3.18. The elements lf a for o G O form a basis ofW. 

3.4 G-models are pseudo-toric - step 4 

Let us define a basis on W e made of vectors lf a . If we consider the inclusion 
map: 

in the basis made respectively of lf o and l x , than this is just copying some 
coordinates, because each lf o = ^ l x and o determines uniquely all x- Namely 
the coordinate of If , is copied to all coordinates l Xe , for \e G o e . One can 
also say that the coordinate of l Xe is equal to the coordinate of I f , where o e 

is the orbit of Xe- Of course this shows that the map from riees^C^e) to 
P(Wl) that parameterizes the model is also given by monomials - these are 
exactly monomials from section 2, where we just make some variables equal 
to each other. We may also look at the following commutative diagram: 

UeeEnw e ) nw E ) -■» nw L ) 
lit 

Uee E nWHe) -+ nW HE ) — » F(W L ) 
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This proves the main theorem of this section: 

Theorem 3.19. Let H be a normal, abelian subgroup of a group G C S n , 
that acts transitively and freely on a set A of n elements. Let W be the 
space of matrices invariant with respect to the action of G and let W be the 
vector space spanned freely by elements of A. Then a G-model of (T, W, W) 
is pseudo-toric for any tree T. 

4 Polytopes of G-models 

In this section we will show how the construction from the previous section 
works on Kimura models and we will present the algorithm for constructing 
a polytope of a model for a given group G with a normal subgroup H. The 
method was described in a different language in |SSj . The main difference 
(apart from the notation) is that the authors assumed the existence of a 
friendly labelling function, that described which characters are identified. In 
case of G-models we know precisely what this function is: it associates to a 
given character its orbit of a G action. It is easy to see that this is a friendly 
labelling. 

If G = H this is particulary easy. The polytope has got |G|'' B ' _ ' JV ' vertices 
and the algorithm works in time 0(|7V| (|(7|'' B ' _ ' iV ')) assuming that we can 
perform group operations in unit time. 

Program 1. 1. Orient the edges of the tree from the root. 

2. For each inner vertex choose one outcomming edge. 

3. Make a bisection b : G — > B C 7) G \ where B is the standard basis of 

4- Consider all possible associations of elements of G with not-chosen 
edges (there are |G ? |' f; ' _ ' JV ' such associations). 

5. For each such associations, make a full association by assigning an 
element of G to each chosen edge in such a way that the (signed) sum 
of elements around each inner vertex gives a neutral element in G. 

6. For each full association output the vertex of the polytope: (b(g e ) e& E), 
where g e is the element of the group associated to edge e. 
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Example 4.1. For a 3-Kimura model on a tree with one inner vertex the 
vertices of P correspond to triples of characters of the group that sum up to 
a neutral character: 



1) (0,0), (0,0), (0,0) 

4) (i,o), (i,o), (o,o) 

7) (0,1), (0,1), (0,0) 
10) (1,1), (1,1), (0,0) 
1^(1,0), (1,1), (0,1) 
16) (1,1), (1,0), (0,1) 

This in the coordinates of the lattice gives us vertices of the polytope: 
1) 1,0,0,0,1,0,0,0,1,0,0,0 2) 1,0,0,0,0,1,0,0,0,1,0,0 



2) (0,0), (1,0), (1,0) 
5) (0,0), (0,1), (0,1) 
8) (0,0), (1,1), (1,1) 
11) (0,1), (1,0), (1,1) 

14) (i,o), (o,i), (i,i) 



3) (1,0), (0,0), (1,0) 
6) (0,1), (0,0), (0,1) 
9) (1,1), (0,0), (1,1) 
12) (0,1), (1,1), (1,0) 
15) (1,1), (0,1), (1,0) 



3) 0,1,0,0,1,0,0,0,0,1,0,0 
5) 1,0,0,0,0,0,1,0,0,0,1,0 
7) 0,0,1,0,0,0,1,0,1,0,0,0 
9) 0,0,0,1,1,0,0,0,0,0,0,1 
11) 0,0,1,0,0,1,0,0,0,0,0,1 
13) 0,1,0,0,0,0,0,1,0,0,1,0 
15) 0,0,0,1,0,0,1,0,0,1,0,0 



4) 0,1,0,0,0,1,0,0,1,0,0,0 

6) 0,0,1,0,1,0,0,0,0,0,1,0 
8) 1,0,0,0,0,0,0,1,0,0,0,1 
10) 0,0,0,1,0,0,0,1,1,0,0,0 
12) 0,0,1,0,0,0,0,1,0,1,0,0 

14) 0,1,0,0,0,0,1,0,0,0,0,1 

16) 0,0,0,1,0,1,0,0,0,0,1,0 



The basis for W for 3-Kimura (in previous notation vectors l x = ^2, x(h a 1 ^'>) 
b) is the following: 



1111 
1111 
1111 
1111 



1-1-11 

-11 1-1 

-11 1-1 

1-1-11 



1-11-1 
-11-11 

1-11-1 
-11-11 



1 -1 

1 -1 

-1 1 

-1 1 



For the 2-Kimura model the elements of H are in the order: 

(1)(2)(3)(4);(1,2)(3,4);(1,3)(2,4);(1,4)(2,3) 

and G is spanned by H and (3,4). 

Now if we consider the action of G on H we obtain: 
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• The orbit of the trivial character is of course only the trivial character. 
This tells us that the first vector is in Wq and will be considered as the 
first basis vector. 

• The orbit of the character that associates —1 to (1, 3) (2, 4) and (1, 4) (2, 3) 
has got also only one element. For example let us notice that 

X ((3,4)(l,3)(2,4)(3,4)) = 4)(2, 3)) = -1 = X ((l, 3)(2, 4)). 

This means that the last vector also will be a basis vector of Wq- 

• There is also one more orbit that contains two left characters. If we take 
their sum (as functions, not characters) we obtain a function that asso- 
ciates 2 to (1)(2)(3)(4), -2 to (1,2)(3,4) and to other two elements. 
This gives us an element: 

"2-20 

-2 2 

2 -2 

0-22 

This is of course the sum of two l x . 

We obtain fi = h, f2 = h + h, fz — U, where U are matrices introduced 
above. Let F = {fx, f'2, fs} and L = {li, . . . , Z 4 }. From the previous section 
F is the basis of Wq and L of Wh- This can be checked directly in this 
example. Let us now look at the map for a tree Y with one inner vertex. 
Elements of Wq are special elements of Wh- We have a map: 

(J/ l )j=l,..,3,i=l,...,3 ->■ (^)j=l,...,4,i=l,..,3- 

Here j parameterizes base vectors and i parameterizes edges. Our model is 
the composition of this map and a model map for H. The image of the first 
map is a subspace given by a condition that the coordinates corresponding 
to 1% and 1% are equal for each i = 1, . . . , 3. Let us see this directly: 
The fixed bijection b is the following: 

6(e) = (1,0,0,0), 6(x 3 ) = (0,l,0,0) 

b( X i) = (0,0,1,0), 6(x 2 ) = (0,0,0,l) 
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where Xi an d X3 are i n the same orbit. Now the domain of ip for the group 
H is {(xi, . . . , Xi 2 ) : Xi G C} in the order corresponding to 14.11 (we fix an 
isomorphism with xi — (1, 0) and X3 — (0, 1)). This tells us that the subspace 
n ee s(^ / c?)e is given by conditions x 2 = x 3 (the coordinates of / 2 and /3 for 

Wr), Xq = X 7 , Xio = x n . 

This procedure works generally. After having fixed the polytope for a 
subgroup H, that is in the lattice M (whose coordinates are indexed by 
edges and characters of H) we consider a morphism form M onto the lattice 
M' (whose coordinates are indexed by edges and orbits of characters of H) 
that just assigns a character to a given orbit. This morphism if of course 
just summing up coordinates that are in the same orbit of the action of G on 
H. The image of the polytope P is a polytope of our model. For 3-Kimura 
we sum up coordinates (order from I4.ip obtaining a polytope for 2-Kimura 
model: 

I) 1,0,0,1,0,0,1,0,0 2) 1,0,0,0,1,0,0,1,0 
3) 0,1,0,1,0,0,0,1,0 4) 0,1,0,0,1,0,1,0,0 
5) 1,0,0,0,1,0,0,1,0 6) 0,1,0,1,0,0,0,1,0 
7) 0,1,0,0,1,0,1,0,0 8) 1,0,0,0,0,1,0,0,1 
9) 0,0,1,1,0,0,0,0,1 10) 0,0,1,0,0,1,1,0,0 

II) 0,1,0,0,1,0,0,0,1 12) 0,1,0,0,0,1,0,1,0 
13) 0,1,0,0,0,1,0,1,0 14) 0,1,0,0,1,0,0,0,1 
15) 0,0,1,0,1,0,0,1,0 16) 0,0,1,0,1,0,0,1,0 
what after removing double entries gives vertices: 
1) 1,0,0,1,0,0,1,0,0 2) 1,0,0,0,1,0,0,1,0 

3) 0,1,0,1,0,0,0,1,0 4) 0,1,0,0,1,0,1,0,0 
5) 1,0,0,0,0,1,0,0,1 6) 0,0,1,1,0,0,0,0,1 
7) 0,0,1,0,0,1,1,0,0 8) 0,1,0,0,1,0,0,0,1 
9) 0,1,0,0,0,1,0,1,0 10) 0,0,1,0,1,0,0,1,0 

5 Normality of G-models 

We have seen that the models associated to a group containing a normal, 
abelian case are pseudo-toric. To obtain results that some varieties are toric 
one needs to prove normality. We will see that in general one cannot expect 
a G-model to be normal, but in many cases it is. The normality of varieties 
is equivalent to normality of polytopes. First let us start with a technical 
lemma. 
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Lemma 5.1. Let Pi and P 2 be two normal polytopes contained respectively 
in lattices L\ and L 2 spanned by the points of the polytopes. Suppose that we 
have got morphisms pi : Li —> L of lattices for i — 1,2 such that Pi(Pi) C S, 
where S is a standard symplex (convex hull of standard basis). Then the fiber 
product Pi x l P 2 is normal in the lattice spanned by its points. 

Proof. Let q G n{P\ x l Pi) and let qi be the projection of q to Lj. One can 
see that as q was the sum of points that belong to Pi Xl P2 with coefficients 
summing up to n and was in the convex hull of n times the points of Pi XlP 2 , 
then each is the sum of points that belong to P; with coefficients summing 
up to n and is in the convex hull of n times the points of Pj. This means 
that qi G nPi. From the assumptions we obtain: 

n 

with each G p. We also know that Pi(qi) = ^2(^2) and this is an element 
of nS. Moreover Pi(vj) G S. But let us notice that each element of nS can 
be uniquely written as the sum of n elements of S. This means that the 
collections (pi(v\) , . . . , pi(vD) and (p 2 (v\), . . . ,p 2 (v^)) are the same up to 
permutation, so we can assume that Pi(vj) = p 2 (v?). So we can lift each pair 
(vj,v?) to a point Vj G Pi x L P 2 that projects respectively to vj and u|. One 
obtains q = Y^j=i v i wna t finishes the proof. □ 

Let us consider two trees Ti and T 2 (in an abelian case) and two leaves 
li G i — 1,2. Let Pi G Mi be a polytope for a tree Tj for i — 1,2. Now 
we can take the lattice L whose basis is given by characters of G. We take 
the projection p^ : Mj — > L onto the coordinates indexed by pairs (li, x) x eG*- 
From the construction of the polytope of the model we know that its vertices 
correspond to networks. The fiber product corresponds to such pairs of 
networks that associate the same character to li and l 2 . This is the same as 
associating a character to each edge of the tree T, that is obtained from Pi 
and T2 by gluing chosen leaves. We see that the fiber product of polytopes 
gives a polytope of a new tree. 

This is also true in a general (not necessarily abelian) case. Now the 
vertices correspond to associations of orbits of characters of H to each edge, 
such that we can find for each edge a representative and the representatives 
form a network. Let us notice that by the conjugation action of G on H* we 
can always choose a representative on one fixed leaf arbitrarily. This means 
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that if the orbits associated to leaves l\ and l 2 agree than, we can find a 
network on the tree T that projects to a fixed element of the orbit. This 
proves that also in a general case the polytope of the tree T is given as a 
fiber product of polytopes for trees Tj. 

Remark 5.2. Of course one can notice that the above assumptions are equiv- 
alent to the fact that the labelling function from |SS] was friendly. 

The projections Pi of the polytopes Pi are indeed in S. This means that 
thanks to lemma 15.11 we only have to consider normality of a model for a 
3-leaf tree, as each trivalent tree can be constructed by gluing such trees. 

Proposition 5.3. The G-models for abelian groups: Z 2 , Z 2 x Z 2 , Z 3 and Z 4 

are normal (for trivalent trees). 

Proof. One can easily find the polytopes for these models with the 3-leaf tree 
using the algorithm from the previous section. One can find the coordinates 
in the lattice that is bigger than the lattice spanned by the points of the 
polytope, but after an easy change of coordinates one can find the right 
presentation of the polytope. We can also compute the Grobner basis of its 
cone and see that it is in the polytope (we use computer program Macaulay) 
what proves normality. □ 

Remark 5.4. The polytope of a 2-Kimura model is not normal. 

Proof. The polytope of a 2-Kimura model is: 
1) 1,0,0,1,0,0,1,0,0 2) 1,0,0,0,1,0,0,1,0 
3) 0,1,0,1,0,0,0,1,0 4) 0,1,0,0,1,0,1,0,0 
5) 1,0,0,0,0,1,0,0,1 6) 0,0,1,1,0,0,0,0,1 
7) 0,0,1,0,0,1,1,0,0 8) 0,1,0,0,1,0,0,0,1 
9) 0,1,0,0,0,1,0,1,0 10) 0,0,1,0,1,0,0,1,0 

The point (1,0,1,1,0,1,1,0,1) is in the lattice spanned by the vertices 
of the polytope and in IP but it is not sum of two points of P. □ 



6 Open Problems 

From the previous section we see that all checked G-models for an abelian 
group G are toric. It would be interesting to characterize groups for which 
G-models are toric. It is natural to ask whether each G-model for an abelian 
group is normal. 
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Appendix 



Here we show an explicit example when the equality of the parameters before 
the Fourier transform does not imply the equality after it. 

Let us consider the group G = Z 4 . The space of matrices is of the form: 

abed 
d a b c 
c d a b 
b c d a 

The matrix of the type above corresponds to a function / : G — > C, such that 
/(0) = a, f(l) = b, f{2) = c and /(3) = d. Now the Fourier transform of / 
gives us: f(xo) = a + b + c+ d, f(xi) = a + ib- c-id, f{x2) = a-b + c- d, 
fiXs) = a — ib — c + id. If we consider a submodel defined by /(l) = 
/(2) what corresponds to b = c the Fourier transform gives us respectively 
(x , Xi, x 2 , x 3 ) = (a + 2b + d, a + (i — 1)6 — id, a — d, a — (i + l)b + id), what 
is a linear subspace defined as (i + l)x\ — 2ix 2 + (i — l)x 3 = 0, hence is not 
an equality of two distinct variables. 
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